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ABSTRACT 

We report the discovery of a candidate z ~ 8 galaxy lensed by a foreground z ^ 1.7 
group, found in the pure-parallel Hubble Space Telescope images taken as part of the 
BoRG Survey. The system consists of a lo98-dropout galaxy with ttiab — 25.92 ±0.14 
which appears to be strongly lensed by a foreground group with a photometric redshift 
of z ~ 1.7, which would be the highest-redshift deflector discovered to date. Several 
pieces of circumstantial evidence point toward the lensing hypothesis: 1) The z ~ 8 
dropout is close in projection to the foreground group, 2) it is slightly elongated tan- 
gentially in the direction of the putative critical curve, 3) a possible counter image 
is marginally detected and 4) the >098-dropout is not clustered, implying that it is 
intrinsically fainter that it appears, i.e. it has been magnified. We show that the con- 
figuration can be accurately reproduced by a Singular Isothermal Ellipsoid deflector 
and a source with a Sersic surface brightness profile. The best fit model has a deflec- 
tor Einstein Radius of l'/49±Hi> 

and the unlensed z ~ 8 source has a Sersic index 

n = 0.92±g;^ and an effective radius of r e = O'/OS^;^ with m AB = 27 A ±0.1. The 
magnification of the source at the location of the dropout is 3. 7 ±0.2. A non-parametric 
reconstruction of the source with LENSVIEW corroborates the derived source parame- 
ters. The lens model yields a group mass of 1.081q!i5 x 10 12 ^© and a mass-to-light 
ratio of M/Lb = 30.7^126^©/-^© within the Einstein Radius. Assuming the stan- 
dard dark matter fraction found for lower redshift lenses, the main deflector galaxy 
has a stellar mass-to-light ratio of M*/Lb = 2.2^' ^Mq/Lq. If confirmed the system 
provides a unique opportunity to obtain spectroscopy of a typical z ~ 8 galaxy and an 
absolute measure of mass for a galaxy at z ~ 1.7. The discovery of this candidate is 
qualitatively consistent with the theoretical expectation that ~10%ofz~8 galaxies 
with luminosity L > 2L* are significantly magnified. However, deeper Hubble obser- 
vations and spectroscopic confirmation are required to rule out alternative scenarios 
and confirm the lensing hypothesis. 

Key words: galaxies: high-redshift, cosmology: observations, gravitational lensing: 
strong 



1 INTRODUCTION 



The epoch of reionisation at redshift z ~ 7 — 12 (Komatsu 
et al.|2010||Shull et al.|2012 ) is a key period in the evolution 



of the Universe when the first galaxies formed and reionised 
the neutral Hydrogen in the intergalactic medium (IGM). 



Due to the absorption of emitted photons with energies 
high enough to ionise neutral Hydrogen, the spectra of these 
early galaxies have highly attenuated regions bluer than the 
Lyman-break at 1216A, enabling their identification from 
broad band imaging with the Lyman-Break, or dropout, 



technique (Steidel et al. 1996 Giavalisco 2002). Until re- 
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cently, this technique had been successfully used to build 
large galaxy samples out to z ~ 6, with sources selected us- 
ing observations from both from ground- and space-based 
observations. However, the sample sizes were dramatically 
reduced at z > 7, because of relative inefficiency of near- 
infrared imaging (Bou wens &; Ilhngworth|2007 Oesch et al. 



2008). The introduction of Wide Field Camera 3 (WFC3) 



on the Hubble Space Telescope (HST) has provided new op- 
portunities in the field, with the construction of the first 
sizeable samples of galaxies at a photometric redshift z > 7, 
with over 200 such objects identified to date (Bu nker et al.| 

201011^ i - -hi— i r^— 



2011a 



McLure et al.||2010| |Yan et al.||2010| Bouwens et al.l 



stein et al.||2012| |Ellis et al.||2012| |Oesch et al.||2013| ). 

These new observations have opened the possibility of 
quantifying and studying the galaxy luminosity function, 
4>(L), during the epoch of reionisation as well as its evolu- 
tion with redshift. Measurement of the luminosity function 
is a key requirement for understanding the production of 
ionising photons at high-z which completes hydrogen reion- 
isation. In this context, it is of fundamental importance to 
discriminate between intrinsic evolution of the galaxy lumi- 
nosity function and observational bias. 

One possible effect is gravitational lensing magnifica- 
tion bias (e.g., Wallington &; Narayan 1993). As suggested 
by |Wyithe et al. (2011), this can be significant given the 
luminosity function at z > 7, even for fields without mas- 



sive structures such as the Hubble Ultradeep Field (Beck- 
with et al.||2007 ). Typically, strong gravitational lensing of 



high-z galaxies is achieved around the critical lines (i.e. high- 
magnification regions) of intermediate redshift massive clus- 
ters 



(e.g., see 



Bradac et al. 2009 Bradley et al. 2012a). 



Through the use of these gravitational telescopes, it is pos- 
sible to achieve magnifications of « 10 and higher (e.g. Yee 
et al. 1996}, thus reaching fainter luminosities and detect- 



ing intrinsically ultrafaint galaxies with high efficiency. This 
is illustrated by a recent report of a candidate at z ~ 12 



strongly lensed with multiple images ( |Coe et aL 2013) as 
well as a handful of other candidate strongly lensed galaxies 
at z ~ 8 and above ( |Zheng et al.|2012" Bouwens et al.|2012| ). 
Gravitational lensing also allows the possibility of spectro- 
scopic follow-up of the candidates ( Treu et al.|2012 ; Bradac 



et al.|[2012 ). In general, magnification bias depends on the 



slope of the luminosity function of the sources: steeper lu- 
minosity functions yield more lensed sources because the 
loss in solid angle is compensated by the steep rise in num- 
ber counts at deeper flux limits. Therefore, this effect be- 
comes especially important for the bright end of the galaxy 
luminosity function (L > L*), where L* is the "knee" of 
the luminosity function in the Schechter parameterization, 
<t>(L) = <t>«(L/L i .) a e X p- (L/L * ) . 

However, magnification is not restricted to massive clus- 
ters. Foreground structures can result in positive magnifica- 



tion along more typical lines of sight. Wyithe et al. (2011) 
showed that ~ 10% of galaxies at z ~ 8 with L > 2L* 
should be magnified by a factor of a few, owing to lens- 
ing by galaxies and groups along the line of sight. To in- 
vestigate the impact of gravitational lensing on the bright- 
end of the observed galaxy luminosity function at z ~ 8, 
we examined the best candidates identified in the Brightest 
of Reionising Galaxies (BoRG) WFC3 HST survey (hav- 
ing vrtAB < 26.5 and S/N > 8). Although in a different 



context, previous studies have already shown that multiply- 
imaged sources are common amongst apparently bright in- 
frared sources found by wide area Herschel surveys ( Ne grello] 
et al.||2010||Gonzalez-Nuevo et aL|2012| ). To investigate the 



impact of gravitational lensing on the bright-end of the ob- 
served galaxy luminosity function at z ~ 8, we examined 
the best candidates identified in the large-area Brightest 
of Reionizing Galaxies (BoRG) W FC3 HST survey ( |Trenti| 
|et al.|2012| |Bradley et al.|2012b| , searching for gravitation- 
ally lensed configurations. 

Throughout this paper we assume a ACDM cosmology 
of On = 0.27, 0\ = 0.73 and H = 70kms" 1 Mpc" 1 . Coor- 
dinates are given in the J2000.0 epoch and AB magnitudes 
are used. This paper is structured as follows. We discuss the 
data, photometry and photometric redshifts in Section [2] 
and the case for lensing in Section [3] In Section [4] we intro- 
duce the Markov-Chain Monte Carlo (MCMC) method that 
we have applied to investigate the lensing hypothesis, as well 
as the MCMC results. In Section [5] we compare the z ~ 8 
source properties to those in the literature. In Section [6] we 
investigate the properties of the deflector, and we conclude 
in Section [71 



SEARCH FOR LENSED Z 
THE BORG SURVEY 



SOURCES IN 



The BoRG survey is a Hubble WFC3 large pure-parallel pro- 
gram aimed at identifying galaxy candidates at z ~ 8 from 



their colours based on the Lyman-break technique (Steidel 
et al.|1996| ). The core of the BoRG survey is based on WFC3 
imaging in four filters (V<506, ^098, J125, Hieo) obtained as 
part of HST programs #11700 and #12572 (PI Trenti), com- 
plemented by archival data with a similar observational de- 
sign and same infrared (IR) filters, but with F600LP substi- 
tuting F606W (HST program # 11702, PI Yan) ( |Yan et al. 
2011| |Yan fe Collaboration||2012) . A fu ll description of the 



survey and of its findings is reported in Trenti et al. (2011 



2012) and |Bradley et al.|p012b) . To summarize, the BoRG 

8 LBGs identified as Y098- 
and 27.4 mag 



survey has so far found 33 z 
dropouts with J125 magnitudes between 25 



(Bradley et al. 2012b). The candidates were obtained by 
analyzing data over a total area of 274 arcmin 2 divided in 
59 independent lines of sight. The large number of inde- 
pendent fields results in a sample unaffected by large-scale 
structure uncertainty and includes a variety of different envi- 
ronments in terms of line-of-sight structure, representing an 
ideal dataset to quantify magnification bias. On average, the 
distribution of magnification along a random line of sight is 
skewed and peaks slightly below one (underdensities) with a 
long tail to positive magnification (overdensities) (e.g., Suyu 



et al.||2010 and references therein). 



To search for evidence of large magnifications, we vi- 
sually inspected the regions surrounding the 8 sources with 
the highest S/N in the survey and identified one candidate 
(borg_0440-5244_682 in the notation of|Bradley et al.|2012b|), 



which is close in projection to a small group of foreground 
galaxies. We also find a possible counter image on the op- 
posite side of the foreground group, making this source a 
good candidate to be magnified significantly and warranting 
further investigation. Should this source be strongly lensed, 
we would find the empirical lensed fraction of the brightest 
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Figure 1. Postage stamp images of the candidate lensed z ~ 8 galaxy in field BoRG_0440-5244 in the four HST/WFC3 filters (V600, 
^098? Jl25, Hiqo from left to right). The images are 4'/ 8 on each side. The location of the Yo98-dropout galaxy is marked by the large 
black circle. The possible faint counter image on the opposite side of the lens is indicated by the small black circle. 



galaxies in BoRG to be ~ 12.5%. 
ment with the predicted value of ^ 



This value is in agree- 
- 10% for a survey such 



as BoRG by Wyithe et al. (2011) 



2.1 



Photometry of Ybs 
sources 



-dropout and foreground 



The candidate z ~ 8 galaxy, borg_0440-5244_682 was iden- 
tified in field borg_0440-5244 (following the notation of 
Brad ley et al.|2012b ). The photometry presented in Bradley 



et al. (2012b) gives a J125 magnitude of 25.9 ±0.1 mag, a 



^098 — ^125 residual of > 2.1 mag, and S/N after sky sub- 
traction in V600, ^098, J125 and Hieo of —0.8, —0.5,9.1 and 
5.7, respectively. We find a J125 ellipticity for the dropout 
of ^ 0.4 in the direction of the suspected critical curve 
of the lens. The J125 half-light radius of the dropout was 
measured, using SExtractor, to be 0V21 (l.lkpc at z ~ 8). 
Within a few arcseconds of the dropout there is a compact 
group of three objects. The two brighter objects have essen- 
tially identical colours and are identified as a single source 
with J125 00 23.5 mag in our SExtractor catalog (source 
borg_0440-5244_647). The third galaxy has J125 ~ 24.6 mag 
(identified as borg_0440-5244_650). The photometry of all 
objects is summarised in Table [I] Postage-stamp images of 
the dropout and foreground galaxies in the single optical 
and three near-IR filters of the BoRG survey are presented 
in Figure [l] Visual inspection of the images highlights a pos- 
sible very faint (S/N~ 2) counter image, which is apparent 
primarily in the deepest band, J125, and located immedi- 
ately below the foreground galaxies opposite the primary 
dropout image. We present the labeled dropout images and 
deflector galaxies in Figure [2] 

To investigate this configuration further, we perform 
aperture photometry in IRAF on the primary and counter 
images to determine the flux ratio in J125 and H\qq. To 
minimize contamination of the counter image photometry 
by the lensing group, we operate on a version of the J125 and 
.Hi 60 images where we subtracted a scaled and PSF-matched 
version of the Y098 data. For this, we use ISIS2.2 ( |Alard 
2000) to produce a convolution kernel between the images 
in the different bands and to match the PSFs. There is some 
residual due to the foreground galaxies having intrinsically 
different profiles between filters. The galaxies on the whole, 
though, are adequately subtracted, and the residuals do not 




COUNTER-IMAGE? 



Figure 2. The lens configuration (J125, 4'/ 8 on a side), with 
labels identifying the lensing group, the primary image of the 
-dropout and the possible counter image. 



affect the photometry at the image positions. The subtracted 
images are shown in Figure [3] 

Photometry is performed in a 4 pixel circular aperture 
(0 ; /32). In J125, the integrated flux in the primary image is 
2678 ± 305 e~ (consistent with the S/N= 9.1 determination 
Brad ley" et al.|[2012b ), and the flux in the counter image 



is 560 ±301 e~. The uncertainty includes both poisson noise 
and background. This gives a flux ratio of 4.8 =b 2.6, where 
the large error propagates from the low signal-to- noise of the 
counter image. 

We do not detect a clear signal at the location of the 
counter image in the iJi60 image, but this is not too surpris- 
ing since from the flux ratio measured in J125 and the H±eo 
primary image 5.7cr detection, we would expect the Hieo 
counter image to be detected only at 1.2 ± 0.7cr. 

In the Markov-Chain Monte Carlo analysis (Section |4|, 
we fit both bands simultaneously so the non-detection will 
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Figure 3. Normalised residuals of the subtractions of the Yo98 image from the J125 image (left) and Hiqq image (right). The images 
are 4 ; /8 on each side. There remains some residual of the deflector galaxy cores due to the varying light profiles of the galaxies between 
filters. We notice the primary image has a clear residual, while we can marginally detect the counter image at a S/N of ~ 2 in the J125 
subtraction. We do not detect the counter image in the Hiqq. The location of the dropout objects are circled in black. See Figure [2] for 
location of individual galaxies/IDs. 



ID 


^600 


^098 


^125 


^160 


682 


> 28.26 


> 28.15 


25.92 ±0.14 


25.75 ±0.19 


647 


24.20 ±0.07 


23.87 ±0.03 


23.41 ±0.02 


23.49 ± 0.04 


650 


25.17 ±0.13 


24.87 ±0.05 


24.58 ±0.14 


24.31 ±0.31 



Table 1. The photometry for the primary dropout (ID: 682), the main lensing group (ID: 647) and a nearby galaxy (ID: 650) which is 
potentially part of the group. 



be taken into account in the modelling and contribute to 
constrain the lens and source properties. 



2.2 Photometric redshifts of foreground galaxies 

To constrain the redshift of the foreground sources, we 
utilise the bayesian photometric redshift code BPZ (|Ben-| 



itez 



2008; |Coe et al.||2007| ). The two brightest galaxies in 



the group (source borg_0440-5244_647) have a redshift dis- 
tribution strongly peaked at z ~ 1.77, primarily determined 
by the assumption that the 4000 A break is the cause for 
the red Y-J colour of the galaxies. The limited information 
from our four-filter photometry does not allow us to con- 
strain accurately the redshift of source borg_0440-5244_650, 
which is characterized by a more extended p(z), peaking at 
z ~ 0.3. However, the probability of this fainter galaxy being 
at z ^ 1.7 is no n- negligible, leading to the likely possibility 
that the broad p(z) is simply an effect of photometric scatter 
and not of a different redshift. To support the group hypoth- 
esis, we searched for galaxies in the immediate proximity of 
the group and dropout sources. In a square field centred on 
the lens candidate with side length 56", we find 10 galax- 
ies brighter than m — 26 mag with photometric redshifts 
peaked at 1.5 < z < 2.1. 

Similar results on the photometric redshift distribution 



of the foreground galaxies are confirmed by the Zebra code 
(Fel dmann et al.| [2006), which prefers a high- z solution for 
borg_0440-5244_650 as well, assuming a continuous star for- 
mation history for that galaxy. If the lens galaxies are all 
at z rsj 1.77, their luminosity translates into an halo mass of 
- 4 x 10 12 M for borg_0440-5244_647 and - 1.5 x 10 12 M 
for borg_0440-5244_650 respectively, on the basis of abun- 
dance matching based on the galaxy luminosity obtained 
using the Cooray (2005) conditional luminosity function. 



At face value this would give circular velocities V c ~ 300 
km s _1 and V c ~ 240 km s _1 respectively for the halos of 
the two galaxies. Translating circular velocities into central 
velocity dispersion ao of a Single Isothermal Sphere model 
for the lens is complicated by baryonic physics which gen- 
erally tend to increase the central velocity dispersion of a 
collapsed Navarro &; White ( 19961) dark matter halo. In prac- 



tice, it is comparable to the observed stellar velocity disper- 



sion cr* w V c at the scale considered here (e.g., |Bundy et al 
2007 Dutton et al.|2010 ), which is in turn very close to the 
velocity dispersion of the Singular Isothermal Sphere that 
describes the total mass density profile in the inner regions 
of massive e llipticals (|Treu &; Koop mans 2004] Koopmans 



et al. 2009 Auger et al. 2010). In addition, lensing opti- 



cal d epth is a strongly rising function of velocity dispersion 



(e.g. Ruff et al. 2011) and ao ~ 300 km s is a typical 
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value for group-sized halos (Treu & Koopmans 2002). In the 
less likely \ow-z solution for the group (z ~ 0.5), the sources 
would only have V c ~ 100 km s _1 , and therefore be too 
small to act as strong lenses. 



3 THE CASE FOR STRONG LENSING OF 

BORG_0440-5244_682 

We see 4 indicators that provide circumstantial evidence 
that borg_0440-5244_682 is a z ~ 8 LBG lensed by a fore- 
ground group. 

(i) We observe a possible faint counter image on 
the opposite side of the deflector. Both the brighter 
image (the 'primary image'), and the fainter image (the 
'counter image') are completely absent in the Yogs and Veoo 
filters, consistent with the z ~ 8 LBG selection criteria, 
and with the lensing hypothesis (lensing is achromatic). 
Multiple imaging by a deflector is a necessary condition for 
strong lensing. 

(ii) A group of massive galaxies at z ~ 1.7 is found 
in between the candidate images. The two dropouts 
are both strikingly close to a foreground group of galaxies. 
The foreground group members have photometric redshifts 
of z ^ 1.7, which is relatively robust despite the availability 
of only four bands due to the Balmer break occurring in 
the imaged Yo98- As we will show below, this configuration 
requires an Einstein Radius of ~ l'/5, which is consistent 
with strong lensing by a massive galaxy/group at moderate 
redshift ( |Treu fe Koopmans||2002| |Treu|2010| ). 



(hi) The dropout is consistent with being elon- 
gated tangentially in the direction of the suspected 
critical curve. From the photometry presented in Sec- 
tion |2.1| we see the primary image shows tangential 
elongation, consistent with the lensing hypothesis. The 
observed elongation is consistent with that predicted by 
the model once smearing by the WFC3 PSF and low S/N 
of the data are taken into account, as described in Section [4] 



(iv) Bright Yo98-dropouts are observed to be 
strongly clustered. We do not observe clustering in 
this case. If the primary observed dropout was not mag- 
nified, it would be intrinsically very luminous. Therefore, 
we we would expect to find additional Yo98 dropouts in the 
field ( |Trenti et aT||2012) . We do not find any oth er y 98 
dropouts in BoRG_0440-5244 ( |Bradley et"aL|2012b) . Lens- 
ing magnification would imply the dropout is intrinsically 
fainter, and hence less strongly clustered, consistent with 
the observed configuration. 



4 MEASUREMENT OF LENS AND SOURCE 
PROPERTIES USING A MARKOV-CHAIN 
MONTE CARLO BAYESIAN INFERENCE 

In this section, we discuss our inference of source and deflec- 
tor properties using the WFC3 data. We have implemented 
the Metropolis algorithm in a Markov-Chain Monte Carlo 



analysis to find the posterior probability distribution func- 
tion of source and deflector parameters within a Bayesian 
framework. Our deflector mass model consists of a Singular 
Isothermal Ellipsoid (SIE). The SIE model is a generaliza- 
tion of the Singular Isothermal Sphere, and is a commonly- 
used parametrisation of the mass distribution of a lensing 
galaxy ( Keeton|2001a ). The SIE is the simplest mass model 
that provides an adequate representation of the mass distri- 
bution of m assive ellipticals (|Kochanek|1995||Treu &; Koop-| 



2004| |Koopmans et al.||2009| |Treu||2010| ). We adopt 



the definitions of the Einstein Radius, ellipticity and posi- 
tion angle as described by Keeton (2001a). The surface mass 



density, in units of critical density, is defined by, 
1 b 



2 7(1 - e)x* + (1 + e )y* 



(1) 



where b is the Einstein Radius, e = and q is the axis 

ratio of the ellipsoid. According to this definition the Ein- 
stein Radius is defined along the intermediate axis of the 
ellipsoid. When calculating the mass enclosed by the criti- 
cal curve, we integrate the surface mass density within the 
Einstein Radius. The Einstein Radius is then related to the 
velocity dispersion of the SIE profile by 

<J v 2D ds 

b = 47T ~ d7' (2) 

where b is the Einstein Radius, a v is the line-of-sight veloc- 
ity dispersion, q is the axis ratio of the ellipsoid, D s is the 
angular diameter distance from the observer to the source 
and Dds is the angular diamter distance from the deflector 
to the source. The centre of the SIE is fixed to coincide with 
the centre of light of the main deflector galaxy. The deflector 
model therefore has three free parameters: Einstein Radius 
(as a proxy for velocity dispersion), ellipticity and position 
angle. 

For the source surface brightness distribution we adopt 
the Sersic profile (Sersic 1963). The Sersic profile is defined 
by six free parameters: the Sersic index, the effective radius, 
flux, x-position, y-position, ellipticity and position angle. 
This gives a total of 10 parameters for the deflector-source 
model. 



4.1 The MCMC method and code 

We wrote an MCMC sampler in python using the Metropolis 
algorithm. The code accepts data in the form of FITS files 
and calculates an image plane model by calling GRAVLENS 
(Kee ton||2001b| ) to produce an image plane surface bright- 
ness map. Other than this dependency, the code relies only 
upon some widely- available python modules. The user of the 
code has full control over the simulation. They may set the 
step size, number of iterations, initial guesses for parameter 
values, which parameters to vary and which to fix at a cer- 
tain value. The code also has the option to self-calibrate the 
step size to achieve a desired success rate. 

The code follows a standard MCMC process. The log- 
likelihood is |x 2 , where X 2 = ( model ;~ data - ) 2 , the data are 
the WFC3 data in Hieo and J125, and the model is the PSF- 
convolved image plane. We use uniform priors on all param- 
eters. As lensing is achromatic, the only difference between 
the Hieo model and the J125 model comes from the PSFs 
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and the noise model. Hence, the log-likelihood is calculated 
for the Hiqo and J125 independently and summed to gener- 
ate a total log-likelihood for a given model. 



4.2 MCMC lens modelling 

We run the MCMC code with the 10 parameters described 
above. The data are cropped WFC3 FITS files in J125 and 
if 160, along with oversampled PSFs and RMS maps. The 
code allows for oversampled PSFs by oversampling the im- 
age plane equivalently. The step size is optimized to deliver a 
success rate of 0.333. As usual, we discard the MCMC burn 
in. We find well constrained posterior probability distribu- 
tion functions (posterior pdf) for each parameter in both the 
source and lens. 

The inferred parameters from the MCMC analysis are 
presented in Table [2] We find a source with a Sersic index 
of n ~ 1 and a major-axis half-light radius of 0.7 kpc (in 
Section [5] we use SExtractor to determine the circularized 
half-light radius for consistency with the literature). The 
apparent magnitude of the source would be 27.4 ±0.1 mag 
if observed unlensed. The magnification at the location of 
the primary image is 3.7 =b 0.2. The total magnification is 
4.4 ±0.2. 

The lens model based on the mean of the posterior of 
each parameter is presented in Figure|4] This model predicts 
a flux ratio of 4.92, which compares well with the observed 
value of 4.8 ± 2.6. The source model is presented in the left 
panel of Figure [5] at a finer pixel scale than the data as 
lensing allows us to recover more spatial information than if 
the source was not lensed. 

We also reconstruct the source in the mean lens model 
using LENSVIEW ( |Wayth fe Webster|2006t . LENSVIEW creates 
a mapping matrix for a given lens model which reverse 
maps image plane pixels to the source plane. We do this 
independently for the J125 and iii60 images. The LENSVIEW 
reconstructions are shown in the centre and right panels of 
Figures [5] and agree well with the MCMC results. 



5 COMPARISON OF THE SOURCE 
MORPHOLOGY WITH THAT OF 
DROPOUTS IN THE LITERATURE 

In this section we compare the predicted source parameters 
with other z ~ 8 galaxies from the literature. The intrinsic 
source profile recovered in our MCMC analysis, and inde- 
pendently using LENSVIEW on the J125 and iii60 data, is 
presented in Figure [5] We use SExtractor on the intrinsic 
unlensed sources to determine their intrinsic effective radii. 
The effective radius definition is not unique in the litera- 
ture. For consistency with previous work on galaxy sizes at 



high-z (Grazian et al. 2012 ) we use SExtractor to determine 
the effective radius, which uses concentric circular aperture 
measurements. In contrast, the effective radius value found 
in the MCMC analysis is defined along the major-axis of the 

Sersic profile, which will vary by a factor of where | 

is the axis ratio. The intrinsic effective radii for the sources 
are presented in Table [3] and is consistent with the results 
of |Oesch et~aLl ( |2009| ), who found an average z ~ 8 LBG 
intrinsic size of ~ 0.5kpc. 

PSF broadening affects the measurement of effective 

radii. To first approximation, r°^ s — ri ntr2 + rp SF , where 

r° bs is the observed effective radius, rl ntr is the intrinsic ef- 
fective radius and rpsF is the full width at half maximum of 
the PSF. Table [3] summarizes the measurements of the ob- 
served size. We compare our results with the size-luminosity 
relation measured at z ~ 7 by |Grazian et al.| ( |2012| in Fig- 
ure [7] after accounting for an expected approximate scaling 
of the galaxy radius with redshift as 1/(1 + z) (Oesch et 
al. 2009; Wyithe & Loeb 2011). We find that the recovered 
properties of our source are consistent with the z ~ 7 sam- 
ple presented in that work. The LENSVIEW reconstruction and 
the source recovered in the MCMC analysis are plotted in 
red in Figure [7] We also plot the size and luminosity of the 
primary image, assuming it is unlensed and find that the 
galaxy would not be an outlier on the relation. Thus, the 
size of the Yo98-dropout is consistent with literature data 
whether it is lensed or not. 



4.3 The image plane 

Based on our lens model we compute the surface brightness 
distribution in the image plane and how it would appear 
in our HST observations owing to the effects of finite 
sampling, resolution and S/N. The results are shown in 
Figure [5] We simulate a noise-less image plane, an image 
plane with the noise and resolution of our data, and the 
image plane with a deeper observation and better sampling. 
They are shown left to right in Figure [6] respectively. 
The simulations show that the current data is consis- 
tent with an elongated image. The right hand panel of 
Figure [6] shows that improving the S/N by a factor of 
two and adopting suitable subpixel dithering to improve 
the sampling will result in a clear detection of elongation 
and of the counter image, if the lensing hypothesis is correct. 



1 Pure-parallel observations, as obtained by BoRG, are not 
dithered because their design cannot affect the primary spectro- 
scopic observations they are associated with. 



6 DEFLECTOR MASS-TO-LIGHT RATIO 

Strong gravitational lensing offers a direct and uniquely pre- 
cise measurement of the deflector total mass distribution. 
Lensing constraints are also free from assumptions of a sys- 
tem's dynamical properties which is critical in other direct 
measurements of a system's mass (|Cappellari et al. 2009). 



Under the assumption that we are, in fact, observing a lens 
system as described in this paper, we can measure the mass 
of a z 1.7 group and its associated most massive galaxies 
with unprecendeted precision. 

To calculate the mass enclosed mass, we integrate the 
SIE mass distribution over the Einstein Radius of 17 49. We 
also fit a de Vaucouleurs profile to the main deflector galaxy 
using GALFIT to determine its effective radius, r e , and its de 
Vaucouleurs luminosity, LdeVauc- The mass within the effec- 
tive radius is also calculated by integrating the SIE mass 
distribution over r e . It should be noted that the effective ra- 
dius of the de Vaucouleurs profile is defined by GALFIT along 
the major axis. For consistency with the mass measurement 
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Parameter 


Posterior median 


68% credible interval 


SIE a 


l'/49 


l'/48 - l'/51 


SIE ellipticity 


0.31 


0.30 - 0.34 


SIE position angle 


-15.3° 


-17.3° - -12.6° 


Source x-position 


cy/39 


0"38 - 0"41 


Source y-position 


0"54 


Of/ 53 - Of/ 5b 


Source flux 


0.20cps 


0.18cps - 0.23cps 


Source ellipticity 


0.70 


0.59 - 0.75 


Source position angle 


-57.0° 


-58.8° - -54.2° 


Source r e 


/ /14 


0'/13 - Of/ 15 


Source Sersic index 


0.92 


0.76 - 1.16 



Table 2. Posterior median and 68% credible interval for the parameters of the deflector and source model. Uniform priors have been 
assumed for all parameters. 




x position (arc seconds) 



Figure 4. The lens model with parameters at the mean of their MCMC distribution. The yellow point represents the source position 
in this model, and the green points are the image positions (for a point source - there is no image plane shape information here). The 
relative size of the green image points represents their magnifications. The red line is the critical curve and the blue lines are the grid 
that GRAVLENS uses. This model produces a flux ratio of 4.92 consistent with the observed value. 



Source 


Angular r* ntr 


Proper r* ntr 


Angular rj 6s 


Proper r° bs 


MCMC analysis 


Of/09 


0.45kpc 


Of/ 16 


0.80kpc 


J125 LENSVIEW 


f /10 


0.50kpc 


Of/ 16 


0.80kpc 


#160 LENSVIEW 


Of/ 10 


0.50kpc 


/ /17 


0.85kpc 



Table 3. The intrinsic and WFC3 observed effective radii of the source recovered in the multi-band MCMC analysis, and reconstructed 
independently for each band using LENSVIEW. 
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0.5" 








Figure 5. Left: The MCMC recovered source, centre: the J125 LENSVIEW reconstruction, right: the H±qq LENSVIEW reconstruction. As 
gravitational lensing magnifies surface area while preserving surface brightness, we can recover spatial information of the source with 
more detail than if it was unlensed. We sample the source plane at a pixel scale 4 times finer than the data. These sources have not been 
convolved with the HST WFC3 PSF. All three panels are 1//2 on each side and shown on the same scale. 




Figure 6. The image plane with varying degrees of noise. On the left is a noise-less image plane at a pixel scale of (y/08/pixel. The 
centre panels hows the image plane including noise with an RMS of the current data at a pixel scale of O'/OS/pixel. The right panels 
shows how the image plane would appear with twice the S/N of the current data drizzled onto a pixel scale of (y/065/pixel. Elongation 
in the primary image is detected, as is the counter image at a high S/N. 



within the Einstein Radius, we integrate the mass over the 
effective radius denned along the intermediate axis. 

Converting the lensing observables and photometry into 
physical masses and luminosities requires knowledge of the 
source and deflector redshifts. Thus, the uncertainty in pho- 
tometric redshifts is the main source of uncertainty for this 
application. In addition, we need to account for system- 
atic errors in photometric measurements. We take all these 
sources of uncertainty into account by performing a Monte 
Carlo simulation to evaluate credible intervals on the de- 
rived parameters. First, we consider the posterior distribu- 
tion for the photometric redshift of the deflector given by 
BPZ, z = 1.77 ±0.27. Second, we set the uncertainty on the 
main galaxy effective radius to be 20%, to account for sys- 
tematic errors which dominate the formal statistical errors 
provided by GALFIT. The error in effective radius propagates 
through the calculation of deflector and main galaxy mass 
as well as the luminosity enclosed within the effective radius. 
Third, we note that the error in the source redshift is not 
as important as the deflector redshift because the angular 
diameter distance is relatively insensitive to redshift around 



z ~ 8. We set this to ±0.5. Finally, we also include errors on 
the solar luminosity, K-correction, systematics in the pho- 
tometry and random errors in the MCMC analysis. The re- 
sulting values for measured quantities are quoted at the 50th 
percentile and the resulting errors on the derived quantities 
are quoted and plotted at the 68% confidence limit. 

The mass enclosed within the Einstein Radius of l ; /49 of 
the SIE mass distribution is 1.08to.i£ x 10 12 M o . Photometry 
of the luminous matter enclosed within the Einstein Radius 
is listed in Table [I] We apply the distance modulus out to 
z = 1.77 and a K-correction using single stellar populations 
models with ages in the range 1-2 Gyr to convert the J125 
magnitude to the rest-frame B band magnitude. We note 
that the redshift B band is very close to the observed J125 
and therefore the K-correction is fairly robust with respect 
to changes in the assumed SED, resulting in a negligible 
additional uncertainty. Combining these measurements, we 
find a total mass-to-light ratio within the Einstein Radius 
of M/Lb — 3O.7ii2 7 6^o/^0 (^ ne uncertainties include the 
effects of the uncertain deflector and sources redshift as de- 
scribed above). This suggests that dark matter dominates 
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Figure 7. Size versus luminosity relation for the Yogs-dropout source under the lensing scenario, as recovered from our MCMC analysis 
and as reconstructed with LENSVIEW (both red). The blue point shows instead the same measurement but assuming that no lensing 
magnification is present. Both scenarios are consistent with the size-luminosity relation measured by |Grazian et ah| ( |2012| ) at z ~ 7, and 
scaled to z ~ 8 assuming r(z) oc 1/(1 + z). The change in luminosity distance from z = 7 to z = 8 is also accounted for. 



the mass of the deflector and supports the argument that 
the deflector is a group. The Einstein Radius is ~ 20 times 
as large as the r e of the main deflector galaxy which is in 
line with a mass distribution dominated by dark matter. 

To obtain the stellar mass-to-light ratio of the central 
deflector, we calculate the lensing mass within r e and ac- 
count for the dark matter fraction as follows. We use GALFIT 
to fit a de Vaucouleurs profile to the main deflector galaxy. 
We find it to have circularised r e = / . / 072±0.014, and a mag- 
nitude of WdeVauc = 23.55 ±0.1 mag. We find the circular 
effective radius to be 0.6212". i^kpc, making the galaxy very 
compact and less luminous than the most massive galax- 
ies at comparable redshift (Da mjanov et alT||20Q9 ) . Within 
the effective radius the SIE mass distribution has a mass 
of 5. 21^4 x 1O 1O M . We again apply the distance modu- 
lus and K-correction to convert the J125 magnitude to the 
rest-frame B band magnitude. We find a mass-to-light ra- 
tio within the effective radius of M/L B = 3.4l° gM /L . 
Adopting the average dark matter fraction of 1/3 found 
for SLACS systems of this mass ( |Treu &; Koopm ans 2004; 



massive galaxies observed at similar redshifts (so-called 'red 
nuggets') flDaddi et al.|2005||Van Dokkum et al.|2008||Dam " 



Auger et al. 2010), we obtain a stellar mass-to- light ratio 
of M+/L = 2.2!^M /L . We find that the central de- 
flecting galaxy is consistent in size and mass with compact, 



janov et al.||2009| |Newman et al.||2012| . The left panel of 
Figure [8] shows the central deflector galaxy on the size-mass 
relation measured by Newman et al. ( |2012| for red nuggets 
at z rsj 1.7. Note however that our mass measurement is 
based on lensing and is therefore insensitive to any system- 
atic uncertainties that may be associated with the inference 
of stellar mass from spectral energy distribution using pop- 
ulation synthesis models. The right panel of Figure [8] com- 
pares our measurement of the stellar mass-to-light ratio to 
previous measurements of the cosmic evolution of the stellar 
mass-to-light ratio obtained with lensing at lower redshifts 
|Treu &; Koopmans] ( |2004| ) and via the evolution of the Fun- 
damental Plane out to z ~ 1.2 |Treu et~aL] ( |2005| ). The re- 
sulting stellar mass-to-light ratio is reasonable and roughly 
consistent with the extrapolation of lower redshift trends. 



7 CONCLUSION 

In this paper we investigate a candidate galaxy-scale gravita- 
tional lensing system with a z ~ 1.7 lens magnifying a^~8 
galaxy. By analysis of the eight brightest z ~ 8 galaxies from 
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Figure 8. Left: Size mass relation at z ~ 1.7 using stellar masses from spectral synthesis models of multicolour photometry ([Newman] 
|et al.|2 012). Our lens candidate (red point) falls right on the relation, using its lensing based estimate of the stellar mass. Right: Cosmic 
evolution of the stellar mass-to-light ratio based on gravitational lenses at lower redshift (Treu & Koopmans 2004), the evolution of the 
Fundamental Plane measured out to z ~ 1.2 ( |Treu et al.|2005| >. The red point represents the candidate lens presented in this paper. 



the HST BoRG survey (reported in [Bradley et al.||2012b] , 
we identified a candidate Yo98-dropout which appears to be 
strongly lensed by a foreground group of massive galaxies 
at a photometric redshift of z ~ 1.7 in field BoRG_0440- 
5244. In Section [3] we argue that the current data can be 
naturally interpreted in the context of strong gravitational 
lensing given that: 

(i) The primary image of the Yo98-dropout (with tuab = 
25.9 ± 0.1) is close in proximity (~ l'/5) to a group of three 
foreground sources with itiab ~ 23.5 — 24.5 that have a 
photometric redshift of z ~ 1.7 and a velocity dispersion of 
~ 300 km/s, as inferred from typical mass- to-light ratios of 
z ^ 1 — 2 galaxies; 

(ii) The primary image is tangentially elongated along the 
proposed critical curve of the deflector; 

(iii) A candidate faint counter image on the opposite side 
of the deflector is detected at a low S/N in J125; 

(iv) Yo98-dropouts with luminosities similar to borg_0440- 
5244_682 have been observed to be preferentially clustered 
( Trenti et al.|2 012), but this does not appear to be the case 
for this field, implying an intrinsically fainter source. 

Quantitatively, we demonstrate in Section [4] that a 
plausible and realistic lens model (SIE deflector and a 
source with Sersic surface brightness profile) reproduces 
all observed data well. The best fit model, derived with 
a MCMC code, predicts a deflector with an Einstein Ra- 
dius of l'.^to oi, and an unlensed z ~ 8 source with Sersic 
index n = 0.92^16 and a circularised effective radius of 
r e = O708±2:oi with m A B = 27 .4 ±0.1. The magnification 
of the source at the location of the dropout is 3.7 ±0.2. The 
lens model yields a group mass of 1.08^o'.i5 x 1O 12 M0 and a 
mass-to-light ratio of M/Lb — 30.7ti 2 7 6 M© / L© within the 



Einstein Radius. Assuming the standard dark matter frac- 
tion found for lower redshift lenses, the main deflector galaxy 
has a stellar mass-to-light ratio of M*/Lb = 2.2l°"^M© / 'L©. 
Overall, the properties of the main lensing galaxy are con- 
sistent with observations of 'red nuggets' at similar redshifts 



flDaddi et al.|2005| |Trujillo et aI.|2006||Buitrago et al.|2008 
|Van Dokkum et al.|2008| ). 

The discovery of the lensed configuration presented in 
this paper is an important contribution to two different ar- 
eas of galaxy formation and evolution. First, it provides the 
first measurement of the fraction of strongly lensed z ~ 8 
galaxies for observations at random pointings (i.e. not clus- 
ters), which we derive to be ~ 12.5%. This ratio is consistent 



with theoretical modelling by Wyithe et al. (2011 ) and high- 



lights that the evolution of the galaxy luminosity function 
at z > 7 might be intrinsically faster than observed at the 
bright-end. Since Wyithe et al. ( 2011| ) showed that the ef- 
fect of magnification bias on the galaxy luminosity function 
increases with redshift, quantifying and verifying its impact 
at z 8 is valuable to construct predictions for upcoming 
and future surveys targeted at finding L > L* galaxies at 
yet higher redshift. Second, we have obtained a lensing mea- 
surement of the stellar mass of a compact galaxy at z ~ 1.7, 
without being affected by any of the systematic uncertain- 
ties that may be associated with the inference of stellar mass 
from spectral energy distribution using population synthesis 
models. This represents a leap forward compared to other 
similar measurements from gravitational lensing, which are 
limited to z < 1 (see Figure |8j). 

The main limitation of the current study is that the 
evidence for the lensing configuration is circumstantial. 
We have clearly shown that the data are consistent with 
the Yo98-dropout being strongly lensed by the foreground 
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sources (likely located at z ~ 1.7) and that this interpreta- 
tion is expected based on modelling of magnification bias 
in high-z large area surveys such as BoRG. However, we 
cannot rule out alternative explanations. The foreground 
group of galaxies could have a lower velocity dispersion 
(for example if it is not at z ~ 1.7 as the photometric 
redshift solution prefers, but rather at z ~ 0.3) and in 
this scenario the tentative detection of the counter image 
would be a noise fluctuation or an unassociated faint 
source. Another alternative to rule out is the possibility 
that the Yo98-dropout primary image is in reality a strong 
emission line source at the redshift of the group, z ~ 1.7, 
with properties similar to the (rare) galaxies observed 
by Atek et al. (2011). However, this latter hypothesis 



would fail to explain the alignment of the Yo98 -dropout 
elongation with the critical curve of the lens and strongly 
contradict the four-band photometric redshift estimates of 
this source. Overall, the lensing configuration derived by 
our work is the most natural interpretation, but clearly 
a validation of the scenario is required. Intriguingly, our 
modelling makes specific predictions testable with deeper 
HST imaging. As shown in Figure [6] improving the S/N 
of the observations by a factor two would both lead to a 
clear detection of the counter image, and to an increased 
ellipticity of the primary image. Such a goal is within 
reach of HST as it would require only a relatively modest 
investment of time (~ lOh of observations). The scenario 
where the Yo98-dropout is an emission line galaxy can 
be ruled out efficiently by spectroscopic follow-up with 
VLT, for example using X-shooter. Such spectroscopic 
follow-up would not only rule out the contamination from 
a z ~ 1.7 emitter, but also provide an opportunity to 
confirm for the first time the redshift of a galaxy at z > 7.5 
by detection of Lya emission. Long slit observations can 
also attempt a concurrent spectroscopic confirmation 
of the redshift of the main deflector while observing the 
Yo98-dropout. This could make galaxy BoRG_0440-5244_682 
not only the most distant spectroscopically confirmed LBG, 
but also the most distant gravitational lens observed to date. 

Acknowledgments This work is partially supported 
by programs number HST/GO 11700 and 12587 from NASA 
through a grant from the space telescope science institute, 
which is operated by the Association of Universities for 
Research in Astronomy, Inc, under NASA contract NAS5- 
26555. JSBW acknowledges the support of the Australian 
Laureate Fellowship. We would like to acknowledge Sharon 
Rapoport for her initial lens modelling of the system. 



REFERENCES 

Alard C., 2000, Astronomy and Astrophysics Supplement 

Series, 144, 363 
Atek H. et al., 2011, The Astrophysical Journal, 743, 121 
Auger M., Treu T., Gavazzi R., Bolton A., Koopmans L., 

Marshall P., 2010, The Astrophysical Journal Letters, 721, 

L163 

Auger M. W., Treu T., Bolton A. S., Gavazzi R., Koopmans 
L. V. E., Marshall P. J., Moustakas L. A., Buries S., 2010, 
ApJ, 724, 511 



Beckwith S. et al., 2007, The Astronomical Journal, 132, 
1729 

Benitez N., 2008, The Astrophysical Journal, 536, 571 
Bouwens R. et al., 2012, arXiv preprint arXiv:1211.2230 
Bouwens R. et al., 2011a, Nature, 469, 504 
Bouwens R. et al., 2011b, The Astrophysical Journal, 737, 
90 

Bouwens R. J., Illingworth G. D., 2007, Highlights of As- 
tronomy, 14, 246 

Bradac M. et al., 2009, The Astrophysical Journal, 706, 
1201 

Bradac M. et al., 2012, ApJL, 755, L7 
Bradley L. et al., 2012a, The Astrophysical Journal, 747, 3 
Bradley L. et al., 2012b, The Astrophysical Journal, 760, 
108 

Buitrago F., Trujillo I., Conselice C. J., Bouwens R. J., 
Dickinson M., Yan H., 2008, The Astrophysical Journal 
Letters, 687, L61 
Bundy K., Treu T., Ellis R. S., 2007, ApJL, 665, L5 
Bunker A. et al., 2010, Monthly Notices of the Royal As- 
tronomical Society, 409, 855 
Cappellari M. et al., 2009, ApJL, 704, L34 
Coe D., Benftez N., Sanchez S., Jee M., Bouwens R., Ford 

H., 2007, The Astronomical Journal, 132, 926 
Coe D. et al., 2013, The Astrophysical Journal, 762, 32 
Cooray A., 2005, Monthly Notices of the Royal Astronom- 
ical Society, 364, 303 
Daddi E. et al., 2005, The Astrophysical Journal, 626, 680 
Damjanov I. et al., 2009, The Astrophysical Journal, 695, 
101 

Dutton A. A., Conroy C., van den Bosch F. C., Prada F., 

More S., 2010, MNRAS, 407, 2 
Ellis R. et al., 2012, The Astrophysical Journal Letters, 

763, L7 

Feldmann R. et al., 2006, Monthly Notices of the Royal 

Astronomical Society, 372, 565 
Finkelstein S. L. et al., 2012, The Astrophysical Journal, 

756, 164 

Giavalisco M., 2002, Annual Review of Astronomy and As- 
trophysics, 40, 579 
Gonzalez-Nuevo J. et al., 2012, ApJ, 749, 65 
Grazian A. et al., 2012, Astronomy & Astrophysics 
Keeton C., 2001a, arXiv preprint astro-ph/0102341 
Keeton C., 2001b, Arxiv preprint astro-ph/0102340 
Kochanek C. S., 1995, ApJ, 445, 559 
Komatsu E. et al., 2010, arXiv preprint arXiv:1001.4538 
Koopmans L. V. E. et al., 2009, ApJL, 703, L51 
McLure R., Dunlop J., Cirasuolo M., Koekemoer A., Sabbi 
E., Stark D., Targett T., Ellis R. S., 2010, Monthly Notices 
of the Royal Astronomical Society, 403, 960 
Navarro J., White S. D., 1996, in SYMPOSIUM- 
INTERNATIONAL ASTRONOMICAL UNION, Vol. 
171, KLUWER ACADEMIC PUBLISHERS GROUP, pp. 
255-258 

Negrello M. et al., 2010, Science, 330, 800 

Newman A. B., Ellis R. S., Bundy K., Treu T., 2012, The 

Astrophysical Journal, 746, 162 
Oesch P. et al., 2009, The Astrophysical Journal Letters, 

709, L21 

Oesch P. et al., 2013, arXiv preprint arXiv:1301.6162 
Oesch P. et al., 2008, The Astrophysical Journal, 690, 1350 
Ruff A., Gavazzi R., Marshall P., Treu T., Auger M., Brault 



© 0000 RAS, MNRAS 000, 000-000 



12 Barone-Nugent et al. 



F., 2011, The Astro-physical Journal, 727, 96 
Sersic J., 1963, Boletin de la Asociacion Argentina de As- 

tronomia La Plata Argentina, 6, 41 
Shull J., Harness A., Trenti M., Smith B., 2012, The As- 

trophysical Journal, 747, 100 
Steidel C, Giavalisco M., Dickinson M., Adelberger K., 

1996, Arxiv preprint astro-ph/9604140 
Suyu S. H., Marshall P. J., Auger M. W., Hilbert S., Bland- 
ford R. D., Koopmans L. V. E., Fassnacht C. D., Treu T., 

2010, ApJ, 711, 201 
Trenti M. et al., 2011, The Astrophysical Journal Letters, 

727, L39 

Trenti M. et al., 2012, The Astrophysical Journal, 746, 55 
Treu T., 2010, ARA&A, 48, 87 

Treu T., Ellis R., Liao T., van Dokkum P., 2005, The As- 
trophysical Journal, 622 

Treu T., Koopmans L., 2002, The Astrophysical Journal, 
575, 87 

Treu T., Koopmans L. V. E., 2004, ApJ, 611, 739 

Treu T., Trenti M., Stiavelli M., Auger M. W., Bradley 
L. D., 2012, ApJ, 747, 27 

Trujillo I. et al., 2006, Monthly Notices of the Royal As- 
tronomical Society: Letters, 373, L36 

Van Dokkum P. G. et al., 2008, The Astrophysical Journal 
Letters, 677, L5 

Wallington S., Narayan R., 1993, ApJ, 403, 517 

Wayth R., Webster R., 2006, Monthly Notices of the Royal 
Astronomical Society, 372, 1187 

Wyithe J., Yan H., Windhorst R., Mao S., 2011, Nature, 
469, 181 

Yan H., Collaboration H., 2012, in American Astronomical 
Society Meeting Abstracts, Vol. 219 

Yan H., et al., 2011, in Bulletin of the American Astronom- 
ical Society, Vol. 43, p. 12805 

Yan H.-J., Windhorst R. A., Hathi N. P., Cohen S. H., Ryan 
R. E., O'Connell R. W., McCarthy P. J., 2010, Research 
in Astronomy and Astrophysics, 10, 867 

Yee H. K. C, Ellingson E., Bechtold J., Carlberg R. C, 
Cuillandre J.-C, 1996, AJ, 111, 1783 

Zheng W. et al., 2012, Nature, 489, 406 



© 0000 RAS, MNRAS 000, 000-000 



